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Abstract 

A uni-directional two-lane road is approximated by a set of two parallel closed 
one-dimensional chains. Two types of cars i.e. slow and fast ones are considered 
in the system. Based on the Nagel-Schreckenberg (Na-Sch) model of traffic flow 
1 15] 5 a set of reaction-diffusion processes are introduced to simulate the behaviour 



of the cars. Fast cars can pass the slow ones using the passing-lane. We write and 
solve the mean field rate equations for the density of slow and fast cars respectively. 
We also investigate the properties of the model through computer simulations and 
obtain the fundamental diagrams . A comparison between our results and v max = 2 
version of Na-Sch model is made. 
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1 Introduction 



In recent years, modeling traffic flow has been the subject of comprehensive studies by 
statistical physicists |]T], [|, |JJ . Needless to say many general phenomena in vehicular 
traffic can be explained in general terms with these models. Distinct traffic states have 
been identified and some of these models have found empirical applications in real traf- 
fic 0, |!| f|. In these investigations, various theoretical approaches namely microscopic 
car-following models || [7j , hydro-dynamical coarse-grained macroscopic models f| || || , 



and gas-kinetic models |T0|, |TT| have been developed in order to find a better quantitative 
as well as qualitative understanding toward vehicular traffic phenomena. Recently as an 
alternative microscopic description, Probabilistic Cellular Automata (PCA) have come 
into play (for an overview see Refs. 



12, 13 



This approach to theoretical description of 
traffic flow is one of the most effective and well-established ones and there is a relatively 
rich amount of results both numeric and analytic in the literature |l|, 14 . 



In PCA models, space (road), time and velocities of vehicles are assumed to take 
discrete values. This realization of traffic flow provides PCA as an ideal tool for the com- 
puter simulation. One of the prototype PCA models is the so-called Nagel-Schreckenberg 



(Na-Sch) model [15] which describes a single-lane traffic flow. Although the initial obser- 
vations of the Na-Sch model were numerical, shortly thereafter analytical technique were 
also proposed |12|, [13|, |l4j]. Analytical treatments to CA are difficult in general. This is 
mainly due to the discreteness and the use of parallel (synchronous) updating procedures 
which produce the largest correlation among the vehicles with regard to other updating 
schemes. 

Soon after its introduction, the Na-Sch model was extended to account for more realis- 
tic situations such as multi-lane traffic flow Jl6, [17|, bi-directional roads |18| and urban 
traffic IJT9, pQ . In multi-lane traffic, fast cars are capable of passing the slow ones by 
using the fast-lane. The possibility of lane-changing allows for these models to exhibit 
non-trivial and interesting properties which are exclusive to multi-lane traffic flow. De- 
spite the quite large approximative methods applied to single-lane Na-Sch based models, 
there are few analytical approaches to multi-lane traffic flow |[21|| . One main reason is the 
large number of rules in PCA modeling multi-lane traffic. In reality, a driver attempting 
to overtake the car ahead (in a uni-directional road) has to take the following criteria 
into consideration: 

1) There must be enough forward space in the passing- lane 

2) There must be enough backward space in passing-lane so that no accident could occur 
between two simultaneously passing cars. 

Moreover, in bi-directional roads additional criteria are necessary for a successful 
passing (for details see [|T^] ). The main purpose of the present paper is to introduce 
an analytical approach to study a uni-directional two-lane road. The approach we use 
is to some extent similar to PCA, however basic differences are distinguishable. The 
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major distinction is concerned with the type of updating scheme. In contrast to PCA 
which are realized in parallel update, our models are based on time-continuous random 
sequential update. The mechanism of modeling the two-lane traffic we use, is based on 
the stochastic reaction-diffusion processes, however the rules have roots in the Na-Sch 
rules. This paper is organized as follows: In section two we define the first model ( model 
I ) and interpret the rules in terms of those in Na-Sch model. Section three starts with 
the Hamiltonian description of the related master equation and continues with mean field 
rate equations and their solutions. The results of the numerical simulation of the model 

1 ends this section. Next we introduce the second model (model II) in section four which 
is formulated in symmetric as well as asymmetric versions and follow the same steps 
performed in section three to obtain the fundamental diagrams of the both versions. The 
paper ends with some concluding remarks in section five. 

2 Definitions of the Models 

In the first model, a uni-directional two-lane road is approximated by a set of two parallel 
one dimensional chains, each with iV sites. The periodic boundary condition applies to 
both. Cars are considered as particles which occupy sites of the chains. Two type of cars 
exist in the system: slow cars which are denoted by A and fast cars denoted by B. Also 
$ represents an empty site. Each site of the chains is either empty, occupied by a slow 
or by a fast car. 

Fast cars can pass the slow ones with certain probabilities while approaching them. The 
bottom lane is the home-lane and cars are only allowed to use the top lane for passing. 
Once the passing process is achieved, they should return to the home-lane. This real- 
ization of a two-lane road is regarded as " asymmetric " type. Nonetheless " symmetric " 
type could also be implemented where passing from the right is allowed as well. 
In model I, we restrict ourselves to " asymmetric " type. The state of the system is 
characterized by two sets of occupation numbers (£1, £ 2) ■ • ■ ; 6v) and (o"i, a 2 , • • • , a^) for 
the home and passing-lane respectively. £j, cr, = 0, 1, 2 where zero refers to an empty site 
whereas one and two refer to a site being occupied by a slow or a fast car respectively. 
To investigate the characteristics of this model, a simplification has been considered. If 
simultaneous two-car occupation of parallel sites of the chains is forbidden, one can de- 
scribe configurations with a single set of occupation numbers {t^} where r» = 0, 1, 2. 



Inspired by the v max = 2 version of the Na-Sch model [jrj, [12[, we propose the follow- 
ing set of stochastic processes which evolve according to a random sequential updating 
scheme: 

A$ -> $A (with rate h) (1) 
-> $5 (with rate p) (2) 
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A$ — > $5 (with rate q) 



(3) 



5$ — > <M (with rate r) 



-> AA (with rate A) 



(4) 
(5) 



BA& -> $.45 (with rate s) (6) 
In order to illustrate the above definitions, let us express their interpretations : 

The first and the second of the above rules correspond to the free moving of slow and 
fast cars respectively. The third one expresses the accelerated movement of slow cars. 
This step corresponds to the so-called acceleration step in the Na-Sch model. The fourth 
rule simulates the behaviour of a driver randomly reducing his/her speed as a result of en- 
vironmental effects, road conditions etc. This step corresponds to the so-called " random 
breaking" step in Na-Sch model. Finally the last two processes simulate the behaviour 
of the fast-car drivers when approaching a slow car. Either they pass the slow car using 
the passing-lane or they prefer to move behind it which give rises to their speed reduction. 

We recall that in Na-Sch model, the forward movement of each car is highly affected 
by the car ahead. Here for simplicity we have considered the two-site interactions and 
only use three-site interaction for the passing process. In this particular case, it is crucial 
that the site ahead of the slow car should be empty. Despite the partial explanation of 
microscopic rules necessary for the description of a traffic flow in a two-lane road, the 
present model ignores the effect of oncoming fast cars (in the passing-lane) on the fast 
car (in the home-lane). In reality, a fast car attempts to overtake provided that there is 
enough back- space behind him in the passing-lane i.e. there is no passing car close to 



him in the passing-lane EBl [T7|. In model (I) passing occurs locally and irrespective of 



the state of passing-lane behind the fast car in the home-lane. 



3 Master Equation and Mean-Field Rate Equations 

The processes (1) to (6) could be regarded as a two-species one- dimensional reaction- 
diffusion stochastic process. This is an example of hard-core driven lattice gas far from 
equilibrium which has proven to be excellent systems for theoretical investigations of 
low dimensional systems out of thermal equilibrium. A large variety of phenomena had 
already been described by driven lattice gases ( for an overview see [^, and the refer- 
ences therein). Using the rates given by (1-6), one can rewrite the corresponding master 
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equation as a Schrodinger-like equation in imaginary time. 



0_ 

di 



\p(t)) 



n\ P (t)) 



(7) 



The explicit form of H could be written down via the rate equations. Let (n k , A ) 
{{n k>B )) denotes the probability that at time t, the site N = k of the chain is occupied by 
a slow (fast) car. The Hamiltonian formulation of master equation allows for evaluating 
the average quantities in a well-established manner. It could be easily verified that the 
following rate equations hold for the average occupation probabilities. 

^( n k,A) = h{n k _ ljA e k ) + r(n k _ liB e k ) + X(n k>B n k+1A ) - h(n k>A e k+1 ) - q{n kjA e k+1 ) (8) 
In the above equation e k stands for 1 — n k}A — n k}B . Similarly for (n ktB ) we have: 

-J t ( n k,B) = g(n fe _i iJ 4efe)+p(n fe _i i Befc) + s(n fc _2,B^fc-i,Aefc)-r(nfc )B e fc +i)-p(n fei Be fc+ i) (9) 

-X(n kjB n k+ljA ) - s(n kjB n k+ljA e k+2 ) 

Apparently the total number of neither slow nor fast cars are conserved according to the 
dynamics and therefore the right hand side of eqs (8,9) cannot be written as a difference 
of two currents. However, the total number of cars i.e. the sum of slow and fast cars is 
a conserved quantity and the time rate of changing (n Atk ) + (n Btk ) is equal to a differ- 
ence between oncoming and outgoing currents. Summing up eqs (8) and (9) yields the 
following discrete form of the continuity equation: 



In which the explicit form of (J k ut ) is given below. 

(J° ut ) = h{n kA e k+1 ) + r{n k:B e k+1 ) + q{n kA e k+1 ) +p(n KB e k+l ) + s{n KB n k+1A e k+2 ) (11) 

Equations (8), (9) and (11) are valid for arbitrary time t, however our particular interest 
is focused in the longtime behaviour of the system where stationarity is established. In 
the steady state regime, one and two-points correlators in (8,9) will be time- independent. 
Equation (10) implies that in steady state the current would be site-independent as ex- 
pected. 

So far, our result have been exact and no approximation has been implemented. At this 
stage and in order to solve equation (8-11) we resort to a mean-field approximation and 
replace the two point correlators with the product of one-point correlators. Moreover, 
since the closed boundary condition has been been applied, it can be anticipated that 
the steady values of (n k , A ) s and (n k ^ B ) s be site-independent and therefore we omit the 
site-dependence subscripts from equations (8-11). Denoting the steady values of (n A ) s 



d 



[<n M > + (n k , B )} = (JD - (Jn 



(10) 
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and (jib)j by n A and n B respectively, the steady current J turns out to be 



J = (hn A + m B + qn A + pn B + sn A n B )(l - n) (12) 
In the above expression, the total density of the cars has been taken to be n 

n A + n B = n (13) 

our final aim is to write J in terms of total density n and the rates . This is performed if 
one writes n A as a function n and the rates. By applying the mean-field approximation to 
the equation (9) in its steady state form and using (13), one obtains the following equation 

r(n - n A )(l — n) + X(n - n A )n A = qn A (l - n) (14) 
which simply yields the solutions: 

1 



nA = 2X 



nX - (1 - n)(q + r) ± (\nX - (1 - n)(q + r)f + 4rn(l - n)Xj : 



(15) 



the solution with the minus sign is unphysical (n A < 0) so the unique solution is the one 
with the positive sign. We Remark that within the mean-filed approach, one also can 
solve the time- dependent version of the equations (8,9). In this case, the equation for 
(n A ) turns out to be: 

^(n A ) = rn{l -n) - [(q + r)(l - n) - nX](n A ) - X((n A )) 2 (16) 
which simply give rises to the following solution: 

n . _ Cie C2 ( C3 -V 

= i -%m < 17 » 

In which C±, C<i and C3 are constants depending on the rates. In the long-time limit, 
the mean concentration of slow cars exponentially relaxes toward the steady value n A . 
Replacing the above n A into the equation (13), one now has the total current J as a 
function of n and the rates. In order to have better insights into the problem, extended 
computer simulation were carried out. Here we present the result of numerical investiga- 
tions of model I. In these computer simulations, the system size is typically 2400. With 
no loss of generality, we re-scale the time so that the rate of hopping a fast car is set to 
one. The speed of slow cars is supposed to be 70 percent of the speed of the fast cars 
which is realized by taking h — 0.7 . The values of q and A are set 1 and 0.7 respectively. 
One sub-update step consists of a random selection of a site, say N = i and developing 
the state of the link + 1) according to the dynamics. One update step contains L 
sub-updates . The typical number of updates developed in order that the system reaches 
to stationarity is 400000 and the averaging has been performed over 500000 updating 
steps. The initial state of the system was prepared randomly i.e. each site is occupied 
with the probability n. figures below show the result of numerical simulations. 
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Total Density n 



Figure 1: current-density diagram for different values of r. s is set to 0.4 . 
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Figure 2: current-density diagrams for different values of s. r is set to 0.2 . 
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Total density n 



Figure 3: density of slow cars versus the total density for s = 0.4 . 
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4 Model II 



4.1 Asymmetric regulation 

The second model we consider, has less resemblance to the Na-Sch model. Here there is 
no specification of fast and slow cars and only one kind of particle exists in the chain, 
nevertheless the distinction between fast and slow cars is realized by their appearance in 
the passing and home-lanes. In this periodic double-chain model the following processes 
occur in a random sequential updating scheme : 



with rate h 



o o 



with rate a 



with rate g 



o o 



• o o • 

— > with rate b 

• • • • 



As depicted, the " asymmetric " regulation has been adopted so that the top-lane can 
only be used for passing. According to the above rules, once a successful passing has 
taken place, the passing car should return to its home-lane unless the next site in the 
home-lane is already occupied. In this circumstance, it can continue to pass the second 
slow car ( multi-passing ). Each site of the double chain takes four different states but 
according to the above dynamics only three of them appears in the course of time. The 
forbidden state is the one in which the passing-lane site is full and its parallel home-lane 
site is empty. Regarding this fact, we characterize the three allowed states by <3>, A and 
B. $ represents the situation where both parallel sites are empty, A represents the case 
of an occupied site in the home-lane and empty parallel site in passing-lane and finally 
B refers to the case of both parallel sites being occupied. 

This notation yields the following reaction-diffusion processes: 



A $ -> $ A (with rate h) 
Ai^$B (with rate a) 



(18) 
(19) 



B $ -> A A (with rate #) (20) 

(with rate b) (21) 

It is worth mentioning that the above model for a two-lane road is simultaneously 
being considered within the approach of Deterministic Cellular Automata (DCA) |25| . 



4.2 Master equation and mean- field approach 

Similar to the steps performed in model I, one can write the following form of discrete 
continuity equation. 

j t [<t»m> + 2K B )] = (Jf) - (Jr> (22) 

in which 

(Jk Ut ) = h(n kA e k+1 ) + b(n kjB n k+1A ) + g{n kjB e k+1 ) + a(n kA n k+1A ) (23) 



The above expression for (J k ) has a clear interpretation in terms of rules (18-21). In 
steady state, the time dependences in the equation disappear and the current will be 
site- independent. Next we apply the mean-field approximation through which all the 
two-point correlators are replaced by the product of one-point correlators. This leads to 
the following equation for J: 

J = hn A (l -n) + b(n - ^-)n A + g(n - ^)(1 - n) + an\ (24) 
where the relation ^ + n B = n has been used. 

In order to obtain J in terms of total density n and the rates, we must write n A as a 
function of n and the rates. This is done by solving the following equation with its left 
hand side set to zero. 

— ua = ^gnsil — n) — 2an A (25) 
The unique physical solution of the above equation is: 

1 



nA = Ta 



{{g\l - n) 2 + 16an(l - n)g)}* - g{\ - n) (26) 



putting (26) in the eq. (24), the current J is now obtained in terms of n and the 
rates. The result of computer simulations are shown in the following set of figures. Here 
the rate b,g and h are chosen to be 1.0, 1.0 and 0.7 respectively while a is varied. We 
recall that "a" measures the tendency of fast cars to pass the slow ones. The simulation 
specifications are the same as those in model I. 
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Figure 5: current-density diagram for different values of passing rate. 
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Figure 6: density of singly occupied sites versus the total density 
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Figure 7: density of doubly-occupied sites versus the total density. 
4.3 Symmetric Regulation 

Here we allow the fast cars to pass rightward as well. In this case, both the top and 
bottom lanes become identical and fast cars can pass the slow ones irrespective of their 
home-lane. In this symmetric two-lane model, each particle hops one site ahead in its 
home-lane provided that the next site is empty. Otherwise it tries to pass the car ahead. 
This attempt is successful if there is an empty site ahead on the opposite lane. The 
following rules illustrates the model definition. 



with rate h 



with rate h 



with rate g 



o o 



o o o • 

— > with rate g 
mm o m 

The astrix symbol indicates that the process in the opposite lane occurs independently 
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of the configuration of the sites filled with astrix. If we denote the state of two parallel 
sites in which the bottom site is empty and the top one is occupied by B, the state of 
simultaneous occupation of parallel sites by C and adopting the notations $ and A as 
the same in the asymmetric version of the model, then it could easily be verified that the 
forms of the discrete continuity equation and the current are as follows: 

j t «o fc ) + (b k ) + 2(c fc » = (J k -i) - (J k ) (27) 

and 

(Jk,k+i) = h((a k e k+1 ) + (a k b k+1 ) + 2(c k e k+1 ) + (c k b k+1 ) + (b k e k+1 ) + (b k a k+1 ) + (c k a k+1 )) 

+ 9 ((hh+i) + (a k a k+1 )) (28) 

Where (a k ), (b k ) and (c k ) refer to the probabilities that at time t, the site N = k of 
the double-chain has one car in bottom lane, one car in top lane and double-occupancy 
in both lanes respectively. In steady state, the system is both time and site independent. 
Denoting the steady values of (ak), (b k ) and (c k ) by a, b and c, one has the relation: 

a + b . , 

— + c = n (29) 

Moreover, the symmetry between the lanes implies that a = b. The steady value a is 
easily found to be obtained from the following equation: 

( g + h)a 2 = hc(l -n) (30) 
Solving the steady-state equation for a, one finds: 

{[h 2 {l - nf + 4/m(l - n)(g + h)])* - h(l - n) 



2(g + h) 

Also equation (31) leads to the following equation for J. 



(31) 



J 



hn(l - n) + h{a 2 + 2a(n - a)} + ga 2 ) (32) 



Where by putting the eq. (31) into it, one reaches to expression for J in terms of n, g and 
h. We remark that the factor two reflects the number of lanes. The result of computer 
simulations are shown in the following set of figures. The value of h is set to one and g 
is varied. 



5 Concluding Remarks 

We have introduced a two-species reaction-diffusion model for description of a uni- 
directional two-lane road. The type of update we have used is random-sequential which 
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Figure 8: current per lane-density diagram for different values of passing rates 
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Figure 9: density of singly occupied sites versus the total density. 
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Figure 10: density of doubly occupied sites versus the total density. 
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Figure 11: number of lane-changing per update versus total density. 
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sounds more appropriate for analytical treatments. In the first model, the result of nu- 
meric simulations are very close to those in mean-field approach which indicates that the 
effects of correlations are suppressed. However in the second model, there are remark- 
able differences between analytical and numeric results, in model I, the current-density 
diagram is slightly affected by changing the passing rate and the passing process has its 
most effect in the intermediate densities. This could be anticipated since in the low and 
high densities the number of passing considerably reduces. The space-time diagrams of 
the model I reveal the discriminating effect of passing. 

In model II (both symmetric and asymmetric), the maximum of J occurs in different 
values of n in simulation and analytical approach. Mean-field predicts a shift toward 
higher densities while in simulation a slight shift toward left is observed. We note that in 
the PCA based models the maximum of J corresponds to a considerable left-shifted value 



of the density [16], [17| . In symmetric version of the model II, we observe an increment of 
the current with regard to the asymmetric version. In contrast to the asymmetric ver- 
sion, the maximum of J in mean-field approach is higher than its value obtained through 
simulation. Although the current diagram (10) appears asymmetrically with respect to 
the density, the lane-changing diagram (13) is symmetric to a high accuracy. 
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